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Abstract 



It is now believed that the limiting distribution function of the largest 
eigenvalue in the three classic random matrix models GOE, GUE and GSE 
describe new universal limit laws for a wide variety of processes arising in 
mathematical physics and interacting particle systems. These distribution 
functions, expressed in terms of a certain Painleve II function, are described 
and their occurences surveyed. 
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1. Random matrix models 

A random matrix model is a probability space (f2,7 ? ,.F) where the sample space 
fi is a set of matrices. There are three classic finite N random matrix models (see, 
e.g. [[3l] and for early history f3^J): 

• Gaussian Orthogonal Ensemble (f3 = 1) 

— f2 = N x N real symmetric matrices 

— V = "unique" measure that is invariant under orthogonal transforma- 
tions and the matrix elements are i.i.d. random variables. Explicitly, the 
density is 



where cat is a normalization constant and dA — Y\ i dAu Yl iK j dAij , the 
product Lebesgue measure on the independent matrix elements. 
• Gaussian Unitary Ensemble ((3 = 2) 
— fl = N x N hermitian matrices 
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— V= "unique" measure that is invariant under unitary transformations 
and the (independent) real and imaginary matrix elements are i.i.d. ran- 
dom variables. 

• Gaussian Symplectic Ensemble (/? = 4) (see for a definition) 

Generally speaking, the interest lies in the -/V — > oo limit of these models. Here 
we concentrate on one aspect of this limit. In all three models the eigenvalues, 
which arc random variables, are real and with probability one they are distinct. If 
Amax(^l) denotes the largest eigenvalue of the random matrix A, then for each of 
the three Gaussian ensembles we introduce the corresponding distribution function 

F N<p (t) := P/3 (A ma x < t) , = 1, 2, 4. 

The basic limit laws @ [47|, [Hj state thatQ 



exist and are given explicitly by 

F 2 (a) = dct(/-A- Ai . 



Ury 



where 



exp ( — y (,t — s)q 2 (x) dx 
Pa{x)M'{y) - Ai'(x)Ai(y) 



ny 



x — y 

acting on L 2 (s, oo) (Airy kernel) 
and g is the unique solution to the Painleve II equation 

q" = sq + 2q 3 

satisfying the condition 

q(s) ~ Ai(s) as s — > oo. 
The orthogonal and symplectic distribution functions are 

Fi(s) = cxp^J^ q(x)dx^ (F 2 ( S )) 1/2 , 

F 4 (s/V2) = cosh (l[ 9{x)dx^ (F 2 (s)) 1/2 . 

Graphs of the densities dFp/ds are in the adjacent figure and some statistics of Fp 
can be found in the table. 



x Here <r is the standard deviation of the Gaussian distribution on the off-diagonal matrix 
elements. For the normalization we've chosen, a = l/\/2; however, for subsequent comparisons, 
the normalization a = y/N is perhaps more natural. 
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Table 1: The mean (/X/3), standard deviation (erg), skewness (Sp) and 
kurtosis (Kp) of Fp. 





M/3 


a l3 






1 


-1.20653 


1.2680 


0.293 


0.165 


2 


-1.77109 


0.9018 


0.224 


0.093 


4 


-2.30688 


0.7195 


0.166 


0.050 



Probability densities 




-4-2 2 

s 



The Airy kernel is an example of an integrable integral operator [ fL9[ and a general 
theory is developed in p9[ . A vertex operator approach to these distributions (and 
many other closely related distribution functions in random matrix theory) was 
initiated by Adler, Shiota and van Moerbeke (See the review article [[HJ for 
further developments of this latter approach.) 

Historically, the discovery of the connection between Painlcvc functions (Pm in 
this case) and Toeplitz/Fredholm determinants appears in work of Wu et al. |33| 
on the spin-spin correlation functions of the two-dimensional Ising model. Painleve 
functions first appear in random matrix theory in Jimbo et al. |pof where they 
prove the Fredholm determinant of the sine kernel is expressible in terms of Py . 
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Gaudin [[L3| (using Mehta's 1 30 then newly invented method of orthogonal polyno- 
mials) was the first to discover the connection between random matrix theory and 
Fredholm determinants. 

1.1. Universality theorems 

A natural question is to ask whether the above limit laws depend upon the un- 
derlying Gaussian assumption on the probability measure. To investigate this for 



unitarily invariant measures (/3 = 2) one replaces in (1.1) 

cxp (-tr(A 2 )) -> exp {-tr(V(A))) . 
Bleher and Its Q choose 

V(A) = gA i -A 2 ,g>0, 

and subsequently a large class of potentials V was analyzed by Deift et al. [Q. 
These analyses require proving new Plancherel-Rotach type formulas for nonclassical 
orthogonal polynomials. The proofs use Ricmann-Hilbcrt methods. It was shown 
that the generic behavior is GUE; and hence, the limit law for the largest eigenvalue 
is i<2. However, by finely tuning the potential new universality classes will emerge 
at the edge of the spectrum. For /3 = 1,4 a universality theorem was proved by 
Stojanovic [Q for the quartic potential. 

In the case of noninvariant measures, Soshnikov ]42] | proved that for real symmet- 
ric Wigner matrices^] (complex hermitian Wigner matrices) the limiting distribution 
of the largest eigenvalue is F\ (respectively, F2). The significance of this result is 
that nongaussian Wigner measures lie outside the "integrable class" (e.g. there are 
no Fredholm determinant representations for the distribution functions) yet the 
limit laws are the same as in the integrable cases. 

2. Appearance of Fp in limit theorems 

In this section we briefly survey the appearances of the limit laws Fp in widely 
differing areas. 

2.1. Combinatorics 

A major breakthrough ocurred with the work of Baik, Deift and Johansson Q 
when they proved that the limiting distribution of the length of the longest increas- 
ing subsequence in a random permutation is F2. Precisely, if £n(&) is the length of 
the longest increasing subsequence in the permutation a G Sn, then 



< s -> F 2 (s) 



#1/6 



2 A symmetric Wigner matrix is a random matrix whose entries on and above the main diag- 
onal are independent and identically distributed random variables with distribution function F. 
Soshnikov assumes F is even and all moments are finite. 
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as N —> oo. Here the probability measure on the permutation group Sn is the 
uniform measure. Further discussion of this result can be found in Johansson's 
contribution to these proceedings p6| . 

Baik and Rains ^ |6) showed by restricting the set of permutations (and these 
restrictions have natural symmetry interpretations) F\ and F± also appear. Even 
the distributions F% and Ff |5(J] arise. By the Robinson-Schensted-Knuth corre- 
spondence, the Baik-Dcift- Johansson result is equivalent to the limiting distribution 
on the number of boxes in the first row of random standard Young tableaux. (The 
measure is the push-forward of the uniform measure on Sn-) These same authors 
conjectured that the limiting distributions of the number of boxes in the second, 
third, etc. rows were the same as the limiting distributions of the next-largest, next- 
next-largest, etc. eigenvalues in GUE. Since these eigenvalue distributions were also 
found in p7fl , they were able to compare the then unpublished numerical work of 
Odlyzko and Rains §33 with the predicted results of random matrix theory. Sub- 
sequently, Baik, Deift and Johansson 0] proved the conjecture for the second row. 
The full conjecture was proved by Okounkov |53| using topological methods and 
by Johansson [E3| and by Borodin, Okounkov and Olshanski |hJ using analyti- 
cal methods. For an interpretation of the Baik-Deift-Johansson result in terms of 
the card game patience sorting, see the very readable review paper by Aldous and 
Diaconis [0. 



2.2. Growth processes 

Growth processes have an extensive history both in the probability literature and 
the physics literature (see, e. g. ||l5| , |29| , Q and references therein), but it was only 
recently that Johansson p2| , |26| proved that the fluctuations about the limiting 
shape in a certain growth model (Corner Growth Model) are F%. Johansson fur- 
ther pointed out that certain symmetry constraints (inspired from the Baik-Rains 
work |^|, |6|), lead to F\ fluctuations. This growth model is in Johansson's con- 
tribution to these proceedings [ p6| where the close analogy to largest eigenvalue 
distributions is explained. 

Subsequently, Baik and Rains and Gravner, Tracy and Widom have shown 
the same distribution functions appearing in closely related lattice growth models. 
Prahofer and Spohn |3^, reinterpreted the work of || in terms of the physicists' 
polynuclear growth model (PNG) thereby clarifying the role of the symmetry pa- 
rameter p. For example, (3 = 2 describes growth from a single droplet where as j3 = 1 
describes growth from a flat substrate. They also related the distributions functions 
Fp to fluctuations of the height function in the KPZ equation ]28|, [29) . (The con- 
nection with the KPZ equation is heuristic.) Thus one expects on physical grounds 
that the fluctuations of any growth process falling into the 1 + 1 KPZ universality 
class will be described by the distribution functions Fp or one of the generalizations 
by Baik and Rains 0. Such a physical conjecture can be tested experimentally; 
and indeed, Timonen and his colleagues have taken up this challenge. Earlier 
Timonen et al. J32| established experimentally that a slow, flameless burning pro- 
cess in a random medium (paper!) is in the 1 + 1 KPZ universality class. This 
sequence of events is a rare instance in which new results in mathematics inspires 
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new experiments in physics. 

In the context of the PNG model, Prahofer and Spohn have given a process 
interpretation, the Airy process, of F2. Further work in this direction can be found 
in Johansson [ p5| . 

There is an extension of the growth model in to growth in a random environ- 
ment. In the following model of interface growth in two dimensions is considered 
by introducing a height function on the sites of a one-dimensional integer lattice 
with the following update rule: the height above the site x increases to the height 
above ar — 1, if the latter height is larger; otherwise the height above x increases by 
one with probability p x . It is assumed that the p x are chosen independently at ran- 
dom with a common distribution function F , and that the intitial state is such that 
the origin is far above the other sites. In the pure regime Gravner- Tracy- Widom 
identify an asymptotic shape and prove that the fluctuations about that shape, nor- 
malized by the square root of the time, are asymptotically normal. This constrasts 
with the quenched version: conditioned on the environment and normalized by the 
cube root of time, the fluctuations almost surely approach the distribution function 
¥2- We mention that these same authors in [[l8f find, under some conditions on 
F at the right edge, a composite regime where now the interface fluctuations are 
governed by the extremal statistics of p x in the annealed case while the fluctuations 
are asymptotically normal in the quenched case. 



2.3. Random tilings 

The Aztec diamond of order n is a tiling by dominoes of the lattice squares 
[m,m+ 1] x [£, I + 1], m, n g Z, that lie inside the region {(x,y) : \x\ + \y\ < n + 1}. 
A domino is a closed 1 x 2 or 2 x 1 rectangle in R 2 with corners in 1? . A typical 
tiling is shown in the accompanying figure. One observes that near the center the 
tiling appears random, called the temperate zone, whereas near the edges the tiling 
is frozen, called the polar zones. It is a result of Jockush, Propp and Shor pl[ 
(see also that as n — > 00 the boundary between the temperate zone and the 
polar zones (appropriately scaled) converges to a circle (Arctic Circle Theorem). 
Johansson Eij proved that the fluctuations about this limiting circle are Fi- 



2.4. Statistics 

Johnstone considers the largest principal component of the covariance matrix 
X l X where X is an n x p data matrix all of whose entries are independent stan- 
dard Gaussian variables and proves that for appropriate centering and scaling, the 
limiting distribution equals F\ in the limit n,p — > 00 with n/p -> 7 £ R + . Sosh- 
nikov 1 43 has removed the Gaussian assumption but requires that n — p = 0(p 1 ^ 3 ). 
Thus we can anticipate applications of the distributions Fp (and particularly F{) 
to the statistical analysis of large data sets. 



2.5. Queuing theory 

Glynn and Whitt (b| consider a series of n single-server queues each with unlim- 
ited waiting space with a first-in and first-out service. Service times are i.i.d. with 
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mean one and variance a with distribution V. The quantity of interest is D(k,n), 
the departure time of customer k (the last customer to be served) from the last 
queue n. For a fixed number of customers, k, they prove that 

D(k, n) — n 

converges in distribution to a certain functional Dk of A:-dimcnsional Brownian 
motion. They show that Dk is independent of the service time distribution V. It 
was shown in || [l6| that Dk is equal in distribution to the largest eigenvalue of a 
k x k GUE random matrix. This fascinating connection has been greatly clarified 
in recent work of O'Connell and Yor |}5| (see also ]3r|). 
From Johansson it follows for V Poisson that 

/ D([xn\,n) - gin \ 



as n — * oo for some explicitly known constants ci and C2 (depending upon x) . 
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2.6. Superconductors 

Vavilov et al. ]52| have conjectured (based upon certain physical assumptions 
supported by numerical work) that the fluctuation of the excitation gap in a metal 
grain or quantum dot induced by the proximity to a superconductor is described 
by F\ for zero magnetic field and by F2 for nonzero magnetic field. They conclude 
their paper with the remark: 

The universality of our prediction should offer ample opportunities for 
experimental observation. 

Acknowledgements: This work was supported by the National Science Foun- 
dation through grants DMS-9802122 and DMS-9732687. 
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